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Abstract 

The Variational Method is applied within the context of Supersymmetric Quantum 
Mechanics to provide information about the energy and eigenfunction of the lowest 
levels of a Hamiltonian. The approach is illustrated by the case of the Morse potential 

"Tj" ' applied to several diatomic molecules and the results are compared with stablished 

]Ci ■ results. 

O 

, I. Introduction 

Q>^ \ Two decades ago Supersymmetric Quantum Mechanics, SQM, was born as a simplest 

case of field theory in the study of SUSY breaking mechanism of higher dimensional 
quantum field theories, developed to unify the four fundamental interactions in nature, 
^. namely the electroweak, strong and gravitational interactions. So far, SQM has been 

r-| ! extensively used to explore different aspects of non-relativistic quantum mechanical 

systems (see, for instance, [0]-|0 and references therein). The new algebraic method is 
specially good to treat analytically solvable potential problems, 0,0. Nevertheless, 
the approximation methods have shown to be of great interest to provide new results 
in potential problems. In this context, one has recently suggested that the variational 
method could be applied within the formalism of SQM, 0-|0|. This approach has 
already been successfully applied to the Hulthen potential, 0. Here the technique 
is applied to the Morse potential, one of the most successful models to represent the 
states of diatomic molecules, 0,0. The results are shown to be in good agreamment 
with the numerical data and results from the other approximation techniques. 

The formalism of variational method can be found in several text books of quan- 
tum mechanics, (see for instance JTOf); it can be used to obtain the approximate 



X 



energy spectra of potentials, in particular, the ground-state. The approach consists 
in using the trial wave-function \E'(r) that depends of a number of parameters and in 
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varying these parameters in the expression for the expectation value of the energ; 



E 



;y 



where H is the Hamiltonian of the system. The variational parameters in \E' are varied 
until the expectation value of the energy is minimum. This value is an upper limit of 
the energy level. Usually this method is applied only to get the ground state energy, 
although it can be applied to get the energy of the excited states as well. 

Thus the central point of the variational method approach is the obtention of the 
trial wavefunction. At this crucial point the superalgebra is used to obtain this func- 
tion. In what follows, we briefly review the main principles of SQM and introduce the 
variational method. The case of angular momentum / = 0, which is exactly solvable, 
is explicitly evaluated from SQM formalism in order to get the hierarchy. We then 
apply the variational scheme to the Morse potential for different values of / and for 
various diatomic molecules and compare the results with known results from other 



numerical methods, ITlll, [O . 



II. Supersymmetric Quantum Mechanics 

In SQM, |1[-01) for A^ = 2 we have two nilpotent operators, Q and Q^, that satisfy 
the algebra 

{Q,Q+} = Hss ;g2 = g+' = 0, (2) 

where Hss is the supersymmetric Hamiltonian. This algebra can be realized as 



e= :- . Q*= : „ (3) 





where A"^ are bosonic operators. With this realization the supersymmetric Hamilto- 
nian Hss is given by 



Hss=r. .-.. = n „- . (4) 





H"^ are called supersymmetric partner Hamiltonians and share the same spectra, 
apart from the nondegenerate ground state. Using the super-algebra a given Hamil- 
tonian can be factorized in terms of the bosonic operators. In ^ = c = 1 units, it is 
given by 

H, = -\^^+V,{r) = AtA-,+E^^^ (5) 

where Eq is the lowest eigenvalue and the function Vi{r) includes the barrier poten- 
tial term. The bosonic operators are defined by 



where the superpotential Wi{r) satisfies the Riccati equation 

W^ - W[ = 2^1 (r) - 2E^^^ . (7) 

The eigenfunction for the lowest state is related to the superpotential W as 

^(^Vr) = Nexp(- r WAf)df). (8) 

Jo 

Now it is possible to construct the supersymmetric partner Hamiltonian, 

H, = A^At + i?« = -^^ + liW^ + W[) + E«. (9) 

If one factorizes H2 in terms of a new pair of bosonic operators, A2 one gets, 

H2 = AtA-, + Ef = -1^ + \{Wi - W',) + Ef (10) 

('2) 
where Eq is the lowest eigenvalue of II2 and W2 satisfy the Riccati equation, 

Wl -W2 = 2V2{r) - 2E^^^ . (11) 

Thus a whole hierarchy of Hamiltonians can be constructed , with simple relations 
connecting the eigenvalues and eigenfunctions of the n- members, |l|, ||13|| 



^n — A^ A^ 


+ Er 




A^- ^ (^"^ 
" V2 Vdr 


+ Wn{r) 

) 




(1) _ A+A+ ,/,("+l) 


E(^) = 4" 


^+1) 


^W(r) = Nexp{- 


rwi{f)df). 





(12) 
(13) 

(14) 
(15) 
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In this work our interest is in the ground-state eigenfunction, to apply the vari- 
ational method. As shown above, in SQM formalism the ground state eigenfunction 
can be determined from the superpotential W{r), which satisfies the Riccati equation, 
eq.(|^). Usually, if the potential is not exactly solvable, it is easier to determine an 
approximation for the superpotential in eq.(|^ than an eigenfunction that satisfies 
the Schroedinger equation for the original Hamiltonian. Thus we make an ansatz 
in the superpotential based in the superalgebra information in order to evaluate the 
trial wavefunction that minimizes the expectation value of the energy, eq.(|l|). This 
formalism is applied to the Morse potential. 



III. The Morse Potential 



For the diatomic system, the three dimensional Morse oscillator can be written 
as, 

Vm = Die-^"^''-''^^ - 26-"^'"-''^^) (16) 

where D is the dissociation energy, rg is the equihbrium internuclear distance and a 
is the range parameter. We rewrite the original Schroedinger equation H'^ = E^^ in 
terms of a new varible y, 

(-^ + ^^^ + A'(e~'^'~'^^ - 2e-^y-y^^)\ ^{y) = e^{y) (17) 

where y = ar and the constants are set hke 

2 2mD Ti^a^ 

ye = are, A = — -y , E = e- — (18) 

and the parameter m is the reduced mass of the molecule. 

Case / = 

In this case the Schroedinger equation is exactly solvable and therefore the hierar- 
chy of Hamiltonians can be constructed. For this case the Schroedinger equation is 
reduced to the form 



( --^ + A2(e-2(s/-s/<=) - 2e-(^-^^)) ) m{y) = e^(y) 



(19) 



The associated Riccati equation is then 

Wl-W[ + e^^^ = VM (20) 

with 

Vi{y) = \^{e-^^y-y^^ - 2e~^y-y^^) (21) 

whose solution is given by 

Wi{y) = -Xe-^y-y^^ + (A - 1/2) (22) 

and the ground-state energy is 

e« = -(A-l/2)2 (23) 

To find the second member of the super-family we solve the equation 

W^ + Wi + e^o'^ =V2{r) (24) 

We then find that the superpartner is given by 

W2{y) = -\e-^y-y^^ + (A - 3/2) (25) 



and the ground-state energy is 

ef = -{\-?>/2f. (26) 

The process of factorization can be continued and the whole hierarchy can be evalu- 
ated. The result is 

W';+i(</) = -\e-''-"> + (A - -^) (27) 

4"«) = -(A - ?^f. 

Case / ^ 

In this case an analytical exact soluction cannot be determined. Nevertheless, based 
in the above arguments on how to obtain the trial wavefunction, we make the following 
ansatz for the superpotential 

WAy) = -\e-^y-y^^ - ^^-^ + C. (28) 

y 

The first term is taken from the one-dimensional Morse superpotential, case of / = 0, 
eg . (|22D , |]16| . The knowledge of the second term comes from the study of three- 
dimensional potentials, (see, for instance, assyntotically linear potential, |T^, and the 



truncated Coulomb potential, [T^). The c-number C will be taken as the variational 



parameter in the trial wavefunction. 

The eigenfunction obtained from eq. (§) is then 

^{y)o,e-^'-''-''' y'^'e-'^y (29) 

Using this expression as a trial wavefunction in the variational method we change 
the parameter C by the variational parameter /x, i.e., 

m^ = ^(y, C = ^) oc e-^^"''"'^' y'+^ e-^y. (30) 

The energy is then obtained by minimisation of the energy expectation value with 
respect to /i. Thus, the equation to be minimised is 

^ _ !^'^,{y)[-ip + \\e-^^y-y^) - 2e-fa-^-)) + '■^]^Mdy ,^^, 

!^^,{yYdy ■ ^''^ 

We have used this expression to minimize the energy expectation value of various 
molecules: H2, HCl, CO and LiH. The explicit values of the energy for n = and 
different values of I are shown, for known values of their respective potential param- 



eters, [19]: D, a, r^ and m. 



Table 1. Energy eigenvalues (in eV) for different values oil for H2 molecule, with 

"e = 0, 7416A and m = 0, 50391amu. Comparison is 



1.9426/1- 



D = 4.7446e1/, 



made with results from ref.[12 



1 


Variational results 


Shifted 1/A^ 
expansion results 


Modified shifted 1/N 
expansion results 


Exact Numerical 





-4.4758 


-4.4749 


-4.4760 


-4.4759 


5 


-4.2563 


-4.2589 


-4.2592 


-4.2589 


10 


-3.7187 


-3.7247 


-3.7252 


-3.7242 


15 


-2.9578 


-2.9663 


-2.9670 


-2.9654 


20 


-2.0735 


-2.0839 


-2.0846 


-2.0826 



Table 2. Energy eigenvalues (in eV) for different values of / for HCl molecule, 
with D = 37255cm-\ a = 1.8677A-\ r^ = 1.2746A and m = 0.9801045amM. Com- 
parison is made with results from ref.|TI]. 
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Variational results 


Shifted 1/A^ 
expansion results 


Modified shifted 1/N 
expansion results 





-4.4360 


-4.4352 


-4.4355 


5 


-4.3971 


-4.3967 


-4.3968 


10 


-4.2940 


-4.2939 


-4.2940 



Table 3. Energy eigenvalues (in eV) for different values of / for CO molecule. 



with D = 90540cm" 



2.2994A- 



1.1283A and m = 6.8606719amM. Com- 



parison is made with results from ref.|ll 
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Variational results 


Shifted l/N 
expansion results 


Modified shifted 1/N 
expansion results 





-11.093 


-11.091 


-11.092 


5 


-11.085 


-11.084 


-11.084 


10 


-11.066 


-11.065 


-11.065 



Table 4. Energy eigenvalues (in eV) for different values of / for LiH molecule. 



with D = 20287cm- 



1.1280A- 



1.5956A and m = 0.8801221amM. Com- 



parison is made with results from ref.|TT]. 
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Variational results 


Shifted l/N 
expansion results 


Modified shifted l/N 
expansion results 





-2.4291 


-2.4278 


-2.4280 


5 


-2.4014 


-2.3999 


-2.4000 


10 


-2.3287 


-2.3261 


-2.3261 



We stress that in fact we are dealing with an effective potential Veff that satisfies 
Riccati equation, i.e., 

Veff{y) = Wl-W[ + E{ji) (32) 

where Wi = Wi{fi) is given by eq.(PB[) and fl is the parameter that minimises the 
energy of eq. (pTI) . It is given by 



v; 



eff 



-Xe-(y-y^) + (-xe-(y-y^) + ^ - i±i ) - i±i + E{fL) (33) 



y 



y 



The plots of both Veff and Vm (plus the potential barrier term) in the same graph 
for / = 5 are shown bellow. The upper curve corresponds to K//- We notice that 
for low energies the effective potential is quite similar to the real Morse potential and 
this is why our results are in such a good agreement with the other results. 




Figure 1: Plot of Vm plus the potential barrier term and K//- 



IV. Conclusions 



The SQM formalism was used to explore the Morse potential. The hierarchy of 
Hamiltonians was evaluated for the case of angular momentum / = 0, for which the 
Morse potential is exactly solvable. 

We have shown a new approach through the application of the variational method 
with the SQM formalism to evaluate the eigenvalues of the Morse oscillator for 
different values of the angular momentum /. The determination of the trial wave 
function motivated by SQM was made by an ansatz in the superpotential. The 
application of the method for diatomic molecules indicated that the values obtained 
are in very good agreement with shifted large- A^ tecnique results, (SLNT), modified 
large-iV tecnique results (SLNT) and numerical data. 

We conclude that the conception of the trial wavefunction through the superpo- 
tential is a simple and good alternative procedure that enables to evaluate the energy 
eigenvalues with a reasoble accuracy. This is not a priori obvious and it suggests that 
it is a suitable method to treat other potentials. 
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